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Abstract. In the 1930s Tarski showed that real quantifier elimination 
was possible, and in 1975 Collins gave a remotely practicable method, 
albeit with doubly-exponential complexity, which was later shown to 
be inherent. We discuss some of the recent major advances in Collins 
method: such as an alternative approach based on passing via the com¬ 
plexes, and advances which come closer to “solving the question asked” 
rather than “solving all problems to do with these polynomials”. 


1 Introduction 

Although methods with better asymptotic complexity are known in theory (e.g. 
[GV88| ), the workhorse of implemented algorithms for real geometric reasoning is 
Cylindrical Algebraic Decomposition. This was introduced in [Col75| to produce 
a remotely practicable (complexity “merely” doubly exponential in the number 
of variables) alternative to Tarski’s original method from 1930 |Tar51| . whose 
complexity could not be bounded by any tower of exponentials. Tarski in fact 
set out to solve the quantifier elimination problem for real algebraic geometry 
(Section |4]): given Qk+iXk+iQk+ 2 Xk +2 ■ ■ - . ,Xn), where Qi € {V,3} and 

is a Boolean combination of relations involving polynomials piixi,... ,Xn), 
find an equivalent '^'(xi,..., a;^), where if' is a Boolean combination of relations 
involving polynomials qi{xi, ..., Xk)- In fact, we cannot solve this in the language 
of algebraic geometry: we need semi-algebraic geometry, allowing > as welQ as 
=. The necessity of > follows from the fact of the example 3y x = {x > 

0) V (a; = 0); its sufficiency is the point of Tarski’s work. 

2 Cylindrical Algebraic Decomposition by Projection and 
Lifting 

|Col75) constructs a samplec0 Cylindrical Algebraic Decomposition (CAD) of R" 
which is sign-invariant for the pi, where these words are defined as follows. 

^ Strictly speaking > is sufficient, but implementations always allow > and A- In fact, 
7 ^ is intrinsic to the regular chains approach discussed in Section 
^ The “sampled” nature is implicit in | Col75l and its successors], but the authors find 
it helpful to be explicit about this. 
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Definition 1 (CAD terminology). Note that throughout we are ordering our 
coordinates/variables, so that Xn is the “last coordinate”. 

decomposition: a partition of R" into cells Ci indexed by n-tuples of natural 
numbers (so R” = IJj Q and i 7 ^ j ^ Ci n Cj = 0); 

(seini-)algebraic: every Ci is defined by a finite set of equalities and inequalities 
of polynomials in the Xi, including expressions of the form 

RootOhifiixi^y)) < X2 < RootOf3(f2(xi,y)) ( 1 ) 

(where RootOf 2 means “the second real root, counting from — 00 ”); 
cylindrical: for all k < n, if tt^ is the projection onto the first k coordinates, 
then, for all i,j, T^kiCi) and 7 rfc(Cj) are either equal or disjoint; 
sampled: for each cell Ci there is an explicit point Si G Ci; 
sign-invariant: for the polynomials in on each cell, every pi is identically 
zero, or everywhere positive, or everywhere negative. 

Collins constructed such a decomposition by a process now known (at least by 
our colleagues) as CAD by Projection and Lifting (for more details see [DavlSj 'l. 
The key property in this approach is the following. 

Definition 2. A polynomial p{xi,... ,Xm) is delineabl^ over a region C C 
R™-i if: 

1. the portion of the real variety of p that lies in the cylinder C x R over C 
consists of the union of the graphs (called sections) of some k>0 continuous 
functions 9i < ■ ■ ■ < 9k from C to R and; 

2. there exist integers mi,..., mu > 1 such that, for every point (oi,..., Om-i) 
in C, the multiplicity of the root 9i{ai,... , 0 ^- 1 ) 0 /p(ai,..., 0 ^- 1 ,a^m), 
considered as a function of Xm alone, is mt (and in particular is constant). 

A set of polynomials is delineable over C if each is delineable and if the sections 
are either identical or disjoint. This is actually equivalent to saying that the 
product is delineable. 

Intuitively, if the {pi} are delineable over C, their graphs neither fold nor cross. 

Let Vn be the set of polynomials in <1>, with coefficients from some effectiv^ 
field A C R. Then Collins algorithm proceeds as follows: 

^ There are various, subtly different, definitions in the literature. This one is from 
|McC99 |. 

The literature often stipulates Q or the algebraic numbers A. The real requirement 
is that we can perform all the polynomial algebra we need over K, and that, given 
expressions a,b £ K, we can decide the trichotomy a < b or a = b ox a > b. Once we 
start adding transcendental functions to our language, the effectivity of K becomes 
a major problem, as we run across the usual indecidability results. This is addressed 
in different ways in [AMWOS] and |Vor89IVor92| . 
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Projection: Given some Vk C K[xi,... ,Xk] construct a set Vk-i C K[xi, ..., 
Xk-i] such that, over each cell of a CAD sign-invariant for Vk-i, the polyno¬ 
mials of Vk are delineable. Though the details depend on the algorithm, the 
key ingredients are leading coefficients (where these vanish some 6i tends to 
infinity), discriminants (where these vanish some 6i ceases to have constant 
multiplicity) and resultants (where these vanish, the 6i from different poly¬ 
nomials intersect). 

Repeat until we have the set of univariate polynomials Vi. 

Base case: Given Pi, isolate the A^i real roots of these polynomials in R^, 
and construct a CAD consisting of the Ni roots, and the A^i -I- 1 intervals 
between them (or to the left/right of them all). The sample points for the 
0-dimensional cells are the roots themselves: for the 1-dimensional intervals 
we choose any convenient point, generally rational and with denominator 
the smallest power of 2 we can find. 

Lifting: Given a CAD Dk-i of sign-invariant for Vk-i, construct a CAD 

Dk of R*^, sign-invariant for Vk - For each cell Ci, this is done by substituting 
the sample point Si into Vk, and doing the equivalent of the base case for 
the resulting univariate system (valid across the whole of C, if the projection 
operator provides delineable projection polynomials). 

Repeat until we have the CAD of R". 

If we suppose that Vn contains m polynomials, of degree (in each variable) 
bounded by d, and coefficient length bounded by I (coefficients bounded by 2*), 
then the time complexity is bounded |Col75[ Theorem 16] by 

( 2 ) 

This analysis is very sensitive to the details of the sub-algorithms involved, and 
a more refined analysis of the base case |Dav85| reduces the complexity (though 
not the actual running time) to 

O . 

This improvement might seem trivial, but in fact implies taking the fourth root 
of the m, d part of the complexity. 

A less sensitive property (and one that reflects the cost of using such a 
decomposition) is the number of cells: for the Collins method this is bounded, 
by an analysis similar to jBDE~*~14] . by 

O (m^"(2d)2'3"y (3) 

As is often the case in mathematics, we get more insight if we solve an apparently 
harder problem. |McC84) did this, demanding that the decompositions Dk, k < n 
be, not just sign-invariant, but 

order-invariant for the polynomials in i.e. on each cell, every pi is identically 
zero, and vanishes to the same order throughout the cell, or everywhere 
positive, or everywhere negative. 
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This actually lets his Vk be much simpler than Collins’, with the cost that 
the lifting procedure might fail if some element p of Vk nullifies (is identically 
zero) over some cell in Dk-i- In this case, McCallum says that Vk was not 
well-oriented, and has to either: 

1. proceed by working around the problem or concluding it not relevant. This is 
only possible in certain cases (e.g. the cell is dimension 0) [Bro05| . Otherwise; 

2. revert to Collins’ projection (or a variant due to [HonQOj h or, 

3. add the partial derivatives of p to Vk and resume the projection process 
from there — an operation that to the best of the authors’ knowledge has 
never been implemented, doubtless because of the complicated backtracking 
involved, and the fact that, whereas we only ought to add this polynomial 
in the nullifying region, the design of Collins’ algorithm and its successors 
assume a global set of polynomials at each level. 

“Randomly”, well-orientedness ought to occur with probability 1, but we have 
a family of “real-world” examples (simplification/branch cuts, see [BBDP07] ! 
where it often fails. The analogy of ^ is given by [McC851 Theorem 6.1.5] as 



( 4 ) 


and a recent improved analysis in BDE+ldl (12)] reduces this to 



( 5 ) 


3 CAD by Regular Chains 

This alternative to the traditional computation scheme of projection and lifting 
was introduced in [CMXY09] . then improved in |CM14a| . The method can be 
described as “going via the complexes”, since the authors hrst construct a cylin¬ 
drical decomposition of C", and then infer a CAD of R". They make use of the 
well developed body of theory around regular systems |Wan00j for the work over 
the complexes, and the algorithms are all implemented in the RegularChains 
Librarjj^for Maple, hence our designation: CAD by Regular Chains. 

We first need the following analogue of Definition (not precisely analogous, 
as Dehnitionallows for non-square-free polynomials and this does not). 

Definition 3. Let K C C be an effective field. Let C be a subset of and 

P C K[xi, ..., Xn-i,Xn] be a finite set of polynomials whose main variable really 
is Xn. We say that P separates above C if for each a G C: 

1. for each p G P, the polynomial \Cx^{p) does not vanish at a; 

2. the polynomials p{a,Xn) G C[a;„], for all p G P, are squarefree and coprime. 

Note that the empty set is trivially separable. 
http://www.regularchains.org 
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We then need an analogue of Definition[^for the case of complex space. We follow 
[CM14a| and describe these (complex) cylindrical decompositions in terms of the 
tree data structure they are stored as. 

Definition 4. We define a cylindrical decomposition o/C", and its associated 
tree, by induetion on n. 

Base Either: There is one set Di, the whole of C and V = {Di}; 

Base Or: there are r non-constant square-free relatively prime polynomials pi 
such that Di is the set of zeros of pi, and D^+i is the complement: {x : 

Pi{x)p2{x) ■ --Prix) ^ 0}; V = {Di ,..., Dj ., D^_|_ I }. 

Base Tree: The root and all the Di as leaves of it. 

Induction: LetV be a cylindrical decomposition o/C"“^. For each Di G D', let 
Vi be a non-negative integer, and Pi = {pip, ■ ■. ,Pi,ri} be a set of polynomials 
which separates over Di. 

Induction Either: r = 0 and we set Dip = Di x C; 

Induction Or: we set Dij = {{a, x) : a G Di Apip{a,x) = 0}; 

A,r+i = ^ia,x) : a G Di A lljPijia,x) o|; 

Then: a cylindrical decomposition of C" is given by 

v = {A,,: 1 < i < A'l; i<j<n + 1}. 

Induction Tree: If T' is the tree associated to D' then the tree associated to D 
is obtained by adding to each leaf Di G T' as children all the Di j such that 
1 < j < f i + 1 • 

Unlike Definition the different roots of a given polynomial are not separated. 
Each cell is the zero set of a system of polynomial equations and inequations, 
where the main variables are all distinct: a triangular system [ALM99| . 

Definition 5. Let F be a set of polynomials in k = K[xi,... ,x„]. A cylindrical 
decomposition D is E-invariant if, for each cell D G D and each fi G F, either 
f vanishes at all points of D or f vanishes at no point of D. 

The trivial decomposition, obtained by taking the “either” branch each time, 
with one cell, is 0-invariant. Given a cylindrical decomposition D which is E- 
invariant, and supposing E = E U {/}, [CM14a| shows how to refine E to a 
cylindrical decomposition E which is E-invariant, hence the “incremental” in 
the title of their paper. The key ingredients in this process are again leading 
coefficients, resultants and discriminants. The paper |CMXY0^ shows, assuming 
that E C R, how to construct from E a cylindrical algebraic decomposition of 
R" which is sign-invariant for E. 

The construction of the cylindrical decomposition can be seen, as pointed 
out in |CM14a) . as an analogue of the projection phase of projection and lifting. 
Indeed, if n is small, it is often the case that the polynomials at level i in the tree 
corresponding to E are those in Vn-i- The fundamental difference is that the Vi 
are global structures: over the whole cylindrical algebraic decomposition of R^ 
we need to isolate all the branches of all of Efe+i, whereas there is a tree structure 
underpinning E and the cylindrical algebraic decomposition, which means that 
“polynomials are not considered when they are blatantly not relevant”. 
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Example: Consider the parabola p := ax^ + bx + c and assume the variable 
ordering x c>- b a. Suppose we were to use projection and lifting. Then the 
first projection identifies the coefficients a, b, c and the discriminant with respect 
to x: b^ — 4ac. Subsequent projection do not identify any further projection 
polynomials for this example. Lifting produces CADs sign-invariant for these 4 
projection polynomials, as well as p itself. 

The regular chains approach would start by building the following tree, rep¬ 
resenting a cylindrical decomposition of C": 



This decomposition was produced to be sign-invariant for p. However, it 
does not insist on sign-invariance for the all the other projection polynomials. 
In particular, it is not sign-invariant for b. The polynomial b is included in the 
projection set because its vanishing can determine delineability, but only when 
the coefficient of the higher degree terms vanish. So, when a = 0 it is important 
to ensure b is sign-invariant, but not otherwise. Hence the tree above is doing 
only what is necessary to make the final conclusion about p. 

The next step is to apply real root isolation, extending this tree to one rep¬ 
resenting a CAD. At the top level the case a ^ 0 must split into the two pos¬ 
sibilities: a < 0 and a > 0. For brevity we display only the branch for a < 0 
below (where ri and r 2 represent the two real roots of p in the case where the 
leading coefficient is negative and the discriminant positive). The full tree has 
27 leaves, thus representing a CAD with 27 cells. This compares with a minimal 
CAD of 115 cells produced by projection and lifting to be sign invariant for all 
projection polynomials. 


a < 0 
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Are there significant savings in general? We refer the reader to [BDE+ldl 
Table 1]. Here PL-CAD refers to our implementation of McCallum’s algorithm 
of Sectio n RC -INC-CAD refers to the algorithm of |CM14aj (Section]^; and 
Qepcad |feo03| is another, highly optimised, implementation of McCallum’s al¬ 
gorithm. Where both terminate, Qepcad and PL-CAD often, though not always, 
have the same cell count. RC-INC-CAD does sometimes have the same count, but 
on other examples such as BC-Phisanbut-4, needs only 2007 cells, while both 
implementations of McCallum’s algorithm need 51,763. 

4 Quantifier Elimination 

The original motivation for |Col75| was the following problem. 

Problem 1 (Quantifier Elimination). Let Qi G {3, V}, and Urcf be the language 
of Boolean-connected equalities and inequalities concerning polynomials in AT [a; i, 
..., Xn \, where K is an effective field with Q C AT c R . Given a statement 
(known as a Tarski statement, or a Tarski sentence if fc = 0) 

^ ■ ■ ■ Qn^nf^is^l •>••••> ^n) • 4^ ^ ■^RCF j (6) 

the Quantifier Elimination problem is that of producing an equivalent 

W := ... ,Xk) ■■ i/; G £rcf- (7) 

In particular, fc = 0 is a decision problem: is true? 

If we have a CAD of R" (noting that the Xi must be ordered in the same 
way in Definition and formula ([^) sign-invariant for the polynomials of 
then constructing E is conceptually easy. 

1. The truth of in a cell of is that of (j) at the sample point Sj. 

2. projects to a CAD of R'=. 

3. The truth of ^ in a cell Dj G is then the appropriate (V for 3 etc.) 
Boolean combination of the truth of (j) in the cells of V that project to Dy 

4. E is then the disjunction of the defining formulae for all the Dj for which E 
is true. 

There is a problem in practice with the last step, first pointed out in [Bro99) . In 
the lifting stage, we produce branches 9i of polynomials, with descriptions such as 
“that branch of p(xi ,..., xi) which, above the sample point s = (oi,..., a/_i), 
has the (unique) root in {(3,^)”, and this is not a statement of Trcf- We could 
equally describe it as “the third real branch of p{xi,... ,xi) above s”, but again 
this statement is not in £rcf- Now by Thom’s Lemma [CR88) . we can describe 
this branch in terms of the signs of p and its derivatives, but, whereas these 
derivatives are in the Collins projection, they are not in the McCallum projec¬ 
tion, or in the tree constructed by the method of Section However, when it 
comes to describing Dy we can just add these (as described in |Bro99| for pro¬ 
jection and lifting and in |CM14b) for regular chains CAD construction). The 
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additional cost is negligible, in particular, we do not need them for projection 
(Section]^, or for tree construction (Section]^. 

Though it may depend non-linearly on polynomial degree etc., this process 
is linear in the number of cells in and produces a disjunction of at most as 
many clauses as there are cells in 


5 Lower Bounds 


This last remark is the basis of the complexity lower bounds in |DH88IBD07| . 
Both constructions use the fact that 


3z 

m VxjTi—iVy 

m —1 


(Z/m—1 — Um ^ ^m—l — ^m) V {ym—1 — ^ ^m—1 — ^m) 

^ Ura—l — 1 ) 


( 8 ) 

encodes ym = Fm-iiFni-iixm))■ Bence applying this construct m — 1 times to 
yi= Fi{xi) gives 

y-m = Fi{Fi{- ■ ■ Fi{{xn)) ■ ■ •)• 

2"*“^ times 

This can then be used to produce expressions with n quantifiers and having 
isolated point solutions, hence needing cells to describe them (the 

0(n) terms are n/3 + 0(1) in [BD07| and n/5 + 0(1) in |DB88| b An example 
which needs 2^°* ' cells for all possible variable orders is also produced in |BD07| . 
along with another which needs 2^ * cells in one order, but a constant number 

in another. Bence the great interest in variable order selection methods for CAD 
IDSS04IEBDW14IBEW+14[ to name a few]. 

The construction in (81 uses both 3 and V in a way that cannot be unnested. 


In fact, it is possible |Gri^| to decide Tarski sentences (i.e. no free variables) with 
a cost that is singly-exponential in n, but doubly-exponential in o, the number of 
alternations of 3 and V in (§. These methods, or any methods singly-exponential 
in n, have, in general, not been implemented, though there has been work on 
the purely existential case (for example |BnnO80 . 


6 Equational Constraints 

The methods described in the previous sections produce decompositions which 
are sign- (or order-)invariant for a set of polynomials. In particular, we can apply 
steps [TJI^ of Section]^ to the same CAD to solve § for any other (j) involving the 
same polynomials. Indeed, as long as the Xi stayed in the same order, we could 
change the Qi as well. |Col98] suggested that we could do better if (j) was of the 
form Pi = 0 A 0', as we would not be interested in the behaviour of polynomials 
in (j)' except when pi = 0. This was implemented in [McC99| . who produced a 
CAD which was sign-invariant for pi, and when pi = 0, sign-invariant for the 
polynomials in (f>'. The main effect of this is to reduce the double exponent n 
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of m in ([^ by 1 , i.e. to take the square root of this term, as shown recently in 

|RDE+14] Vi 4b 

It is worth seeing how this works. Consider 

^ ■■= ifi = 0) A ((/2 > 0) V (/3 > 0)). (9) 

Then a |McC84| -stvle projection ignoring the fact that there is an equation 
constraint would contair0three disc(/i) and three res(/i, fj). However, [McC99) 
observes that we are not interested in / 2 , /a except when fi — 0 , and hence we 
need only consider disc(/i) and res(/i, / 2 ), res(/i, /a), half as many polynomials. 
Consider now 


h ■= iigi = 0) A {g 2 > 0)) V ((ga > 0) A (34 = 0)). 


( 10 ) 


A |McC84) -stvle projection ignoring the fact that there is an equation constraint 
would contain four discriminants and six resultants. Although (10) does not 
contain an overt equational constraint, 4>i ^ {gi = 0) V (34 = 0), which is 
4>i ('?i5'4 = 0)) s-iid so the equational constraint gig^ = 0 is implicit. If we 

study 5134 = 0 A in the style of ([^, and drop trivial resultants, we consider 
disc( 5 i 54 ), res(5i54,52) and res( 5 i 54 , 53 )- Using the multiplicative properties of 
resultants and discriminants (which we would certainly do in practice!), this is 
disc( 5 i), disc( 54 ) and all the resultants except res( 52 ;fl 3 )! he. two discriminants 
and five resultants. 

Intuitively res( 5 i, 53 ) and res( 52 > 54 ) are re dundant, but how do we achieve 
this in general? This was solved in |BDE~*~13] . where, rather than producing a 
sign-invariant CAD, we compute truth table invariant (a TTICAD) for the two 
propositions (51 = 0) A (52 > 0) and (53 > 0) A (54 = 0), i.e. on each cell, each of 
these two propositions is either identically true, or identically false. This process 
does indeed remove these two resultants, so we have two discriminants and three 
resultants. 


Example: Consider (10) with 


5 i := -k 5 ^ - 4, g 2 := (x - - {y + 3), 


and 

53 := (x - 3)2 + (5 - 2), 54 := (x - 6 )^ + - 4. 

Figurej^shows the two dimensional cells produced for both a sign-invariant CAD 
and a truth-table invariant CAD, built under ordering x ^ y. The sign-invariant 
CAD has 231 cells (72 full-dimensional but the splitting of the final cylinder is 
out of view) and the TTICAD 67 (22 full-dimensional). 

By comparing the figures we see two types of differences. First, the CAD of 
the real line is split into fewer cells (there are not as many cylinders in R^). This 
is the effect of the reduction in projection polynomials identified, (less univariate 


It would also have some leading coefficients etc., but these are not the main drivers 
of the complexity in McCallum’s projection. 
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Fig. 1. The left is a sign-invariant CAD, and the right a TTICAD, for (101 with the 
polynomials from the example. 



polynomials with real roots to isolate). The second difference is that the full¬ 
dimensional cylinders are no longer split over the dashed lines. This came from an 
improvement in the lifting phase (discussed in detail in |BDE~*~14 ). It leverages 
the projection theory to conclude that we usually only need to lift with respect 
to equational constraints themselves. 


More recently |BDE~*~14) truth-table invariance has been achieved even when 
there is no implicit equational constraint, as with an example of the form 


((hi = 0)A(h2>0))V(h3>0). 


( 11 ) 


The savings that can be achieved depend on the number of equational constraints 
involved in sub-clauses of the parent formula. 

It is also possible to apply equational constraints in the regular chains tech¬ 
nology view of CAD |CM14aj . again even when there is no global equational 
constraint, as in (111 [BCD^14] . 


7 How Reliable is this? 

Cylindrical algebraic decomposition can be used as tool in program verification, 
as in the MetiTarski tool |Paul2| . This leads to the question: who will verify the 
CAD, or at least the inferences we draw from it? We note that a positive answer 
to a purely existential question (equally, a negative answer to a purely universal 
question) is easily verified since we have a witness. The converse questions are 
essentially questions of refutation, see |JdM12j . Questions involving a mixture 
of quantifiers are much harder. 

Almost all current implementations of CAD are based on computer algebra 
systems, which are generally unverified. We can at least compare, on a fairly 
level playing field, the implementations in Maple of four algorithms: see Table 
The classification of the amount of mathematics involved is subjective, but 
we note that |McC84| . and hence [BDEWl^ . relies on [Zar65IZar75| to ju stify 
the smaller projection set compared with |Col75] . [CM14a| and [BCD~*~14] rely 
on, inter alia, |ALM99] . 

There are two challenges involved in verifying a CAD algorithm. 
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Table 1. Comparison of algorithms 


Algorithm Implementation 

Code Lines 
(above Maple) 

Specialist 

Mathematics 

Col75l 

EWBD14I 

2600 

some 

McC84l 

EWBD14I 

2500 

a lot 

CM14al 

CM14al 

5000 

medium 

BDEW13I 

EWBD14I 

3000 

a lot 

BCD+141 

BCD+14 

5500 

medium 


1. There is a “program verification” question of ensuring that the algorithms 
produce the result that they say they will, i.e. that resultants, discriminants, 
real roots etc. are computed correctly. This is non-trivial, to say the least, sit¬ 
ting on top of an unverified computer algebra system, but should be feasible 
for an implementation based on a sound kernel, such as Coq or Isabelle. 

2. There is a “mathematics verification” question whether the resulting de¬ 
composition is truly sign-/order-/truth table-invariant for the inputs. This 
is where the column labelled “Mathematics” in Table comes in. The only 
attempt to produce verified CADs known to the authors, |CM121 in Coq], is 
based, not on |Col75| and its successors, but rather on |BPR061 chapter 2], 
itself essentially that of |Tar51j . 

2a. There is an interesting tension here between “precomputed” and ad hoc ver¬ 
ification. An implementation based in |McC84j would essentially have to ver¬ 
ify the relevant theorems from |Zar65IZar75) , but these could be imported as 
pre-verified lemmas. An implementation based on |CM14a| would verify that 
in this case we had an appropriate cylindrical decomposition of C" which 
in this case translated to an appropriate cylindrical algebraic decomposition 
of R”. 

8 Final thoughts 

The topics we focussed on in this paper are implemented in Maple: 

— CAD by Regular Chains is implemented in the RegularChains Library. 
A version of this ships with the core Maple distribution while the latest 
version is freely available from http://www.regularchains.org/. 

— The authors’ own work (equational constraints, truth-table invariance, sub¬ 
decompositions) is freely available in a Maple package ProjectionCAD. 
The latest version is available from: http://opus.bath.ac.uk/43911/. 

Other implementations of cylindrical algebraic decomposition include: 

— Mathematica [StrOB] : The commands CylindricalDecomposition and 
Reduce can make use of an underlying CAD implementation. These com¬ 
mands can be exceptionally fast but it can be hard to judge the CAD 
components individually as they are just one of several underlying meth¬ 
ods available and the output is in the form of formulae rather than cells. 
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— Qepcad [BroOd] : a dedicated interactive command-line program available 
from http;//www.usna.edu/CS/qepcadweb/B/QEPCAD.html. One notable 
feature is the SLFQ program which can simplify large quantifier free for¬ 
mulae giving more readable output. Sage now has a Qepcad interface. 

— Redlog |SS03| : this Reduce package implements CAD along with other 
quantifier elimination methods such as virtual substitution. 

— SyNRAC |IYAY13j : a Maple package notable for its symbolic-numeric ap¬ 
proach. An older version is available for free download from: 

http://jp.fujitsu.com/group/labs/en/techinfo/freeware/synrac/ 
with more recent advances part of the wider Todai Robot project. 

The only reported experiments to cover all of these implementations were de¬ 
tailed in Section 4 of |BCD~*~1^ . 

Of course, this paper surveyed only a few of the recent advances in cylindrical 
algebraic decomposition. Others include (but are not limited to): 


The use of certified numerics in the lifting phase to minimise the amount of 
symbolic computation required |Str06IIYAY13| . 

Local projection schemes [Strl4] . generic projection schemes |SS03] and sin¬ 
gle CAD cells |Bro13I.TdM12] . 

Problem formulation for CAD |DSS04IBDEW 131 WEDBT^ (projection and 
liftintrl |EBDW14lEBC+ m (regular chains). These all develop heuristics to 
help with choices, while HEW~*~14] applies machine learning in the form of 
support vector machines to pick a heuristic. 

Work on cylindrical algebraic sub-decompositions, which return only a subset 
of the cells in a full CAD [Sei06| . In |WBDEi4] algorithms are given to 
return cells that lie on a prescribed variety, or have a designated dimension, 
while in |WDEB13] these techniques are combined to solve a motion planning 
problem. Note that if restricting to cells of full dimension then sample points 
can always be chosen to be rational, greatly reducing running time. 


There are numerous unsolved problems, both theoretical and practical. Three 
that stand out to the authors are the following. 

1. There is no complexity analysis of the Regular Chains method (though 
clearly it is subject to the lower bounds in Section]^. 

2. There has been much progress in the last forty years, but implementations (at 
least for systems with alternations of quantifiers) are still doubly-exponential 
in the number of variables while the theory suggests we can do better. 

3. Cylindricity is needed in stepof quantifier elimination, as 3 translates into 

V and V into /\. However, in fact we only need this at the points where 3 and 

V alternate, so we can weaken the definition of cylindricity from being true 
for all TTfe to merely being true for those k where Xk and Xk+i Etre governed 
by different quantifiers (or where Xk is unquantified but Xk+i is quantified, 
a concept we can call block-cylindrical. Unfortunately, we currently know of 
no way of computing a block-cylindrical algebraic decomposition without 
computing the full cylindrical algebraic decomposition first. 
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